Jamming proteins with slipknots and their free energy landscape 



o 
o 

(N 
o 

Q 

o 

m 



Joanna I. Sulkowska^, Piotr Sulkowski^''^''' and Jose N. Onuchic^ 
^ Center for Theoretical Biological Physics, 
University of California San Diego, 
Oilman Drive 9500, La Jolla 92037, 
^ Physikahsches Institute and Bethe Center for Theoretical Physics, 
Umversitat Bonn, Nussallee 12, 53115 Bonn, Germany 



California Institute of Technology, Pasadena, 
* Institute for Nuclear Studies, 
Hoza 69, 00-681 Warsaw, Poland 



CA 92215, 



Theoretical studies of stretching proteins with slipknots reveal a surprising growth of their un- 
folding times when the stretching force crosses an intermediate threshold. This behavior arises as 
a consequence of the existence of alternative unfolding routes that are dominant at different force 
ranges. Responsible for longer unfolding times at higher forces is the existence of an intermediate, 
metastable configuration where the slipknot is jammed. Simulations are performed with a coarsed 
grained model with further quantification using a refined description of the geometry of the slip- 
knots. The simulation data is used to determine the free energy landscape (FEL) of the protein, 
which supports recent analytical predictions. 

PACS numbers: 87.15.ap, 87.14.E-, 87.15. La, 82.37.Gk, 87.10.-l-e 
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The large increase in determining new protein struc- 
tures has led to the discovery of several proteins with 
complicated topology. This new fact has arised the ques- 
tion if their energy landscape and the folding mechanism 
is similar to typical proteins. One class of such proteins 
includes knotted proteins which comprise around 1% of 
all structures deposited in the PDB database P, 0] ■ A 
related class of proteins contains more subtle geometric 
configurations called slipknots d, Q . Recent theoretical 
studies using structure-based models (where native con- 
tacts are dominant) suggest that slipknot-like conforma- 
tions act like intermediates during the folding of knotted 
proteins This entire new mechanism is consistent 
with energy landscape theory (FEL) and the funnel con- 
cept 0, 9- It was shown that the slipknot formation 
reduces the topological barrier. Complementing regular 
folding studies, additional information about the land- 
scape was obtained by mechanical manipulation of the 
knotted protein with atomic force microsco py [9| both 
experimentall y in [l3,[lH and theoretically in 12, M\- 
For example, [I4I it has been showen that unfolding pro- 
ceeds via a series of jumps between various metastable 
conformations, a mechanism opposite to the smooth un- 
folding in knotted homopolymers. 

Motivated by these early results, we now propose a uni- 
fied picture for the mechanical unfolding of proteins with 
slipknots. In this Letter this question is addressed by 
explaining the role of topological barriers along their me- 
chanical unfolding pathways. Supported by our previous 
results that knotted proteins can still have a minimally 
frustrated funnel-like energy landscape, structure-based 
theoretical coarse-grained models are used 15| to ana- 
lyze the behavior of a slipknot protein under stretching. 
Studies are performed for the a/ (3 class protein thymi- 



dine kinase (PDB code: le2i [ITi]). 
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FIG. 1: Dependence of the unfolding times r on the stretch- 
ing force F for le2i (solid line, in red). In this Letter we 
describe this mechanism as a superposition of two unfolding 
pathways: I for small forces (dashed (lower) line, in blue), 
and II for intermidiate and large forces (dashed-dotted (up- 
per) line, green). 



Most of our analysis is based on stretching simulations 
under constant force [3]. The crucial signature for this 
process is the overall unfolding time from the beginning 
of the stretching until the protein fully unfolds. Normally 
one expects that the transition between the native and 
the unfolded basins to be limited by overcoming the free 
energy barrier, which gets effectively reduced upon an 
application of a stretching force. The rate by which this 
barrier is reduced depends on the distance between the 
unfolded basin and the top of the barrier measured along 
the stretching coordinate x. This idea was first devel- 
oped in the phenomenological model of Bell , which 
states that the unfolding time t decreases exponentially 

Fx 

with applied stretching force F as t{F) — tqc '-'b^. A 



refined analysis performed in ref. [19| revealed that this 
dependence is more complicated but still monotonically 
decreasing. 

The unfolding times for le2i measured in our simula- 
tions are shown as the red curve in Fig. [TJ In contrast to 
the above expectations, increasing the force in the range 
3-3. 5e/ A surprisingly results in a larger stability of the 
protein, e is the typical effective energy of tertiary na- 
tive contacts that is consistent with the value e/A~71 
pN derived in [l^ . A solution for this paradox is accom- 
plished by realizing that unfolding is dominated by two 
distinct, alternative routes that are dominant at different 
force regimes. A routing switch occurs when threshold is 
crossed between weak and intermediate forces. At higher 
forces, mechanical unfolding is dominated by a route that 
involves a jammed slipknot. This jamming gives rise to 
the unexpected dependence of unfolding time on applied 
force. Characterizing this mechanism is the central goal 
of this Letter. 




FIG. 2: A slipknot (left) consists of a threaded loop {k\ — k2, 
in red) which is partialy threaded through a knotting loop 
{k2 — fca, in blue). An example of a protein configuration with 
a tightened slipknot is shown in the right panel. 



To describe the evolution of a slipknot quantitatively 
requires a refined description. A slipknot is character- 
ized by the three points shown in Fig. [2l The first 
point ki is determined by eliminating amino acids con- 
secutively from one terminus until the knot configura- 
tion is reached (which can be detected e.g. by applying 
the KMT algorithm [1^). The two additional points, 
^2 and ^3, correspond to the ends of this knot. In the 
native state the protein le2i contains a slipknot with 
= 10, fcs = 128, fca = 298. These three points divide 
the slipknot into two loops, which are called the knotting 
loop and the threaded loop. The former one is the loop of 
the trefoil knot and the latter one is threaded through the 
knotting loop. Unfolding of the slipknot upon stretch- 
ing depends on the relative shrinking velocity of these 
two loops (see Fig. [3]). When the threaded loop shrinks 
faster than the knotting loop, the slipknot unties. In the 
opposite case the slipknot gets (temporarily) tightened 
or jammed, resulting in a metastable state associated 
to a local minimum in the protein's FEL. Upon further 
stretching, this configuration eventually also unties. The 
evolution of both loops of the slipknot is encoded in the 



time dependence of the points fci, fc2, fca, see Fig. |3l 
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FIG. 3: The behavior of the slipknot during stretching (top) 
is determined by the relative behavior of its two loops, en- 
coded in the time dependence of k\, k2 and fca (bottom). If the 
threaded loop shrinks faster than the knotting loop, ki merges 
with k2 (bottom left) and the slipknot untightens (pathway I, 
top left). If the knotting loop shrinks faster, k2 approaches fcs 
(bottom right, ~ 14000r) and the slipknot gets temporarily 
tightened (pathway II, top right). This is a metastable state 
which can eventually untie further stretching , with ki finally 
merging with k2 (bottom right, ~ 19000r). Kinetic stud- 
ies were performed slightly above folding temperature using 
overdamped Langevin dynamics with typical folding times of 
lOOOOr. 



Before discussing the stretching of le2i, we explain why 
a slipknot formed by a uniformly elastic polymer should 
smoothly unfold under stretching. To simplify the discus- 
sion we approximate the threaded and knotting loops by 
circles of size Rt and Rk- These two loops shrink during 
stretching and, when the threaded one eventually van- 
ishes, the slipknot gets untied. If both loops have similar 
sizes, the slipknot is very unstable and unties immedi- 
ately. When the threaded loop is much larger than the 
knotting one, Rt >> Rk, untightening can be explained 
as follows. The elastic energy associated to local bend- 
ing is proportional to the square of the curvature. If the 
loop is approximated by a circle of radius R, then its local 

uals R^^. The total elastic 
From the assumption 
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curvature is constant and eq 
energy is § dsR~'^ ~ R~^ 
Rt » Rk we conclude that upon stretching it is ener- 
getically favorable to decrease Rt rather than Rk- This 
happens until both radii become equal and then, just as 
above, the slipknot gets very unstable and untightens. In 
this discussion we have not yet taken into account that 
when a slipknot is stretched some parts of a chain slide 
along each other. This effect could be incorporated by in- 
cluding the friction generated by the sliding [2^ . But in 
the slipknot the sliding region associated with the knot- 
ting loop is much longer than the region associated to 
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the threaded loop. Thus this effect results in a faster 
tightening of the threaded rather than the knotting loop, 
facilitating even more the untightening of the slipknot. 

The above argument should apply to slipknots in 
biomolecules because they are characterized by a per- 
sistence length that in principle is simply related to their 
elasticity 23] . For DNA this effect is described by worm- 
like-chain models (WLC) [ISl and it has been confirmed 
experimentally. Although WLC models are too simple 
to describe the protein general behavior, they are use- 
ful in some limited applications. Thus at first sight one 
might expect that slipknots in proteins should smoothly 
untie upon stretching. Proteins, however, are much more 
complicated than DNA or uniformly elastic polymers. 
The presence of stabilizing native tertiary contacts leads 
to a jumping character during stretching T^. In addi- 
tion their bending energy is not uniform along the chain 
due to the heterogeneity of the amino-acid sequence. As 
a consequence it turns out that the intuition obtained 
through the above analysis of polymers or WLC models 
is misleading. 
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FIG. 4: Dependence on the applied stretching force of the 
probability of choosing pathway I rather than II (see Fig. [Sj. 
This varying probability leads to the complicated dependence 
of the total unfolding time on the stretching force observed in 
Fig. [11 



Our analysis of the evolution of the cndpoints fci , fc2, ^3 
(Fig. [21 bottom) reveals that for various stretching forces 
unfolding proceeds along two distinct pathways (Fig. [3l 
top). In pathway I the slipknot smoothly unties, which 
is observed for relatively weak forces. At intermediate 
forces pathway II starts to dominate and the knotting 
loop can shrink tightly before the threaded one vanishes. 
In this regime the protein gets temporarily jammed (Fig. 
[21 right), leading to much longer unfolding times (catch 
pathway) . The probability of choosing pathway I at dif- 
ferent forces is shown in Fig. [31 This pathway competi- 
tion explains the nontrivial total unfolding time depen- 
dence observed in Fig. [Tl 

The two different pathways I and II arise from com- 
pletely different unfolding mechanisms. Pathway I starts 
and continues mostly from the C-terminal side, along 
16a, 15/3, 14q;, 13/?, 12(helices bundle), 11a (here the 
number denotes a consecutive secondary structure as 



counted from N-terminal, and a or (3 specifies whether 
this is a helix or a /?-sheet; for more details about the 
structure of le2i see the PDB). This is followed by unfold- 
ing of helices 11a, 10a that allows breaking of the con- 
tacts inside the /3-sheet created by the N-terminal, with 
unfolding proceeding also from the N-terminal. Pathway 
II also starts from the C-terminal but rapidly (as soon 
as helix 15 is unfolded) switches to the N-terminal. In 
this case, differently from pathway I, the /3-sheet from 
the N-tcrminal unfolds even before 13/3. These scenarios 
indicate that the pathway I should be dominant at weak 
forces since they are not sufficient to break the /3-sheet 
during first steps of unfolding. The jammed pathway is 
typical only if stretching forces are sufficiently strong for 
unfolding to proceed from the two terminals of the pro- 
tein. 



A similar phenomenon was firstly proposed in ref. 25 1 
and referred to as catch-bonds. Experimental evidence 
suggesting this mechanism was first observed for adhe- 
sion complexes 2^ 2^- Using AFM, at large forces the 
ligand-receptor pair becomes entangled and therefore ex- 
pands the unfolding time. A theoretical description of 

this mechanism was given in ref. [iilli,!!^. 

The kinetic data can also be used to determine the as- 
sociated free energy landscape (FEL) ^. In an initial 
simplification we associate the barriers along the stretch- 
ing coordinate as the the kinetic bottlenecks during the 
mechanical unfolding event. Generalizing Bell's model, 
a recent description of two-state mechanical unfolding in 
the presence of a single transition barrier has been devel- 
oped in [l^, with the rate equation 



-ag) ' 



(1) 

where ly encodes the shape of the barrier. Here denotes 
the distance between the barrier and the unfolded basin 
(in a first approximation it can be regarded as F indepen- 
dent) and lies on the reaction coordinate along the AFM 
pulling direction. It can be experimentally determined 
by measuring how the stretching force modulates the un- 
folding times T. The height of the barrier is denoted by 
AG. Fig. [1] (unfolding times are given by solid red line) 
shows that this single barrier theory is not sufficient for 
the full range of forces. As described before, in the higher 
force regime, additional basins have to be included in the 
energy landscape. Models with several metastable basins 
have been called multi-state FEL models 3l|. Evidence 



supporting the need of multi-states FEL was confirmed 
by AFM experiments in different systems 33, 2^- 

To construct a multi-state FEL that incorporates two 
unfolding pathways I and II we use a linear combina- 
tion of eq. ([IJ-like expressions with different shapes and 
barrier heights. Each one of them essentially accounts 
for the distinct barrier along a relevant unfolding route. 
Fitting the stretching data to eq. H]) with a cusp-like 
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FIG. 5: Pathway II with two barriers. Left: dependence of the 
unfolding time on the applied force with the data and the fit 
to the formula ([T]) for the first maximum (lower, in green) and 
for the second maximum (upper, in blue). Right: schematic 
free energy landscape for this pathway, with jammed slipknot 
in a minimum between two barriers. 



other proteins with shpknots, e.g. Ip6x. Our results 
provide testable predictions that can now be verified by 
AFM stretching experiments. 
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V = \/2 approximation (another possibility = 2/3 for 
the cubic potential in general leads to similar results [Toj ) 
determines accurately the location and the height of the 
potential barriers. Pathway II involves two barriers: first 
until the moment of creation of the intermediate which 
is followed the untieing event. They are characterized by 
(xijAGi) and (a;2,AG'2) arising respectively from the 
lower and upper fits in Fig. [5] (left). The superposition 
of these two fits gives the overall mean unfolding time for 
pathway II (dotted-dashed curve in green in Fig. [1]). For 
the ordinary slipknot unfolding (pathway I), the results 
xi and Gi arise from the dashed blue curve in Fig. [T] 
This analysis leads to the results 



xi = 2.3- 



kuTk 



X2 - 0.7- 



kuTk 



XI 



1.4- 



.kuTk 



AGi = SMbT, AGs = 4.2fcBr, AG/ = 4.7fcBT. 

We conclude that the free energy landscape consists of 
two "valleys" . The force-dependent probability of choos- 
ing one of the valleys during stretching depends on the 
details of the protein structure. It is determined from our 
simulations as shown in Fig. 31 Using these probability 
values and the parameters above for x and AG, we can 
accurately represent the simulation data using a linear 
combination of equations of the form ([T]). This agreement 
supports our analytical analysis and generalizes eq. ([T]) 
for the full of range forces. In addition it demonstrates 
that structure-based models sufhciently capture the ma- 
jor geometrical properties of a slipknotted protein. A 
schematic representation of the free energy landscape for 
pathway II is shown in Fig. [S] (right). 

Summarizing, we have analyzed the process of tighten- 
ing of the slipknot in protein le2i and determined the cor- 
responding free energy landscape. Its main feature is the 
presence of a metastable configuration with a tightened 
slipknot, which is observed for sufficiently large pulling 
forces. This phenomenon does not exist for uniformly 
elastic polymers. In this Letter we concentrated on pro- 
tein le2i but similar behavior has also been observed for 
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